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Abstract. Until now X-radiation from nine millisecond 
pulsars has been detected. In the present paper we sum- 
marize the observations and show the results of a re- 
analysis of archival ROSAT data. In addition we present 
the results of recent observations of PSR J0437— 4715 
with the ROSAT PSPC and HRI detectors. We show 
that the pulsed fraction is independent of energy in the 
range 0.1—2.4 keV. The pulse width as measured at X- 
ray energies is comparable with that observed in the 
radio domain. An upper limit for the X-ray luminos- 
ity of the pulsar's bow-shock visible in Ha is found to 
he Lx < 2 X 10^^ erg s~^(d/18Qpc)^. We further report 
on the discoveries of PSR J1024-0719, J1744-1134 and 
J2124-3358 with the ROSAT-HRI, making them the first 
solitary galactic millisecond pulsars detected at X-ray 
energies. The pulse profile of PSR J2124— 3358 shows 
marginal evidence for a double peak structure. We fi- 
nally discuss the observed emission properties of the de- 
tected millisecond pulsars and conclude that the measured 
power-law spectra and pulse profiles together with the 
close correlation between the pulsar's spin-down energy E 
and the observed X-ray luminosity suggests a non-thermal 
origin for the bulk of the observed X-rays. 

Key words: Pulsars: individual (PSR B1821-24, 
J0218+4232, B1957-K20, J0437-4715, J1012H-5307, J1024- 
0719, J1744-1134, J2124-3358, J0751-I-1807) - X-rays: gen- 
eral - Stars: neutron - Stars: binaries: general - Clusters: 
globular 



1. Introduction 

The satellite observatories ROSAT, ASCA and EUVE 
have brought important progress in neutron star and pul- 
sar astronomy. With significantly higher sensitivities com- 
pared with previous X-ray satellites they allowed for the 
first time the detection of X-ray emission from objects as 
faint as millisecond pulsars. 



Millisecond pulsars form a separate group among the 
rotation-powered pulsars. They are distinguished by their 
small spin periods (P < 20 ms) and high spin stability 
(dP/dt « 10^^^ — 10^^^ s/s). Consequently, they are very 
old objects with spin-down ages P/2P of typically 10^ 

10^" years and magnetic dipol components B± oc 
of the order of 10^ - 10^° G (see [Taylor et al."] l995). More 
than ~ 75% of the known disk millisecond pulsars are in 
binaries with a compact companion star, compared with 
the = 1% of binary pulsars found in the general popula- 
tion. This gives support to the idea that their fast rotation 
has been acquired by angul ar momentum transfer during a 
past mass accretion phase ( Bisnovatyi-Kogan & Romberg 
1974; [Alpar et al. 1982| ; [Bhattacharya fc Van den Heuvcf 
1991; Urpin et al. 1998; see also the recent results by 



Wijnands 



Van der Klis| (1998) on SAX 1808.4-3658). 
launch of ROSAT, nothing was known 
properties of millisecond pul- 



Before the 
on the X-ray emission 
sars. According to the standard models of cooling neu- 
tron stars, millisecond pulsars are too o ld to expect de- 
tectable thermal emission (see Tsuruta 1998 and refer- 
ences therein). Thermal X-ray emission, however, may be 
emitted from polar caps heated up to temperatures of 
few million degrees by high energy secondary stream- 



ing back to the neutron star's polar cap regions (Arons 
1981; iKundt fc Schaa| 1993; |Gil fc Krawc^ 1996). Non- 
thermal processes, which may account for the X-ray emis- 
sion from millisecond pulsars, are: magnetospheric emis- 
sion from relativistic particles, characterized by a power- 
law spectrum; X-rays from a relativistic pulsar wind or 
from an interaction of that wind with the interstellar 
medium or a close companion star ( Arons fc Tavan^ 1993). 

At present, there are 33 rotation- powered pulsars de- 
tected in the soft X-ray domain (see Becker fc Triimpei 
1997 and references therein)^]. Nine of them belong to the 
group of millisecond pulsars. In this paper we review the 
recent millisecond pulsar observations at X-ray energies 
and report on results obtained with the ROSAT, EUVE, 
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ASCA and RXTE satellites. The structure of the paper is 
as follows: After a brief description of the analysis tech- 
niques in Sect. |^ we summarize the X-ray emission prop- 
erties of the faint detections (PSR 1957+20, J075H-1807, 
J1012-F5307, J1024-0719 and J1744-1134) in Sect. |. Re- 
sults from the nearby and bright pulsar PSR J0437— 4715 
are presented in Sect. ^. In chapter ^ and ^ we summarize 
recent results obtained from the globular cluster pulsar 
PSR 1821-24 and from PSR J0218-h4232. The discovery 
of pulsed X-ray emission from the solitary galactic mil- 
lisecond pulsar PSR J2124— 3358 and the conclusions are 
presented in Sect. and Sect. ||. 

2. Data analysis 

The ROSAT data presented in this paper were analyzed 
using the extended scientific analysis software EXSAS 
(Zimmermann et al. 1994). Source positions and count 



pulsation. For instance, evidence for pulsed X-ray emis- 
sion was found in the ROSAT all-sky survey data of PSR 
J0437— 4715 already with a number of ~ 50 source counts 



rates were obtained from a maximum likelihood analy- 
sis of the source photons in combination with a spline fit 
to the background after removing solar-scattered X-rays 
and the particle background. Count rates were dead-time 
and vignetting corrected. 

Five of the nine millisecond pulsars detected by 
ROSAT are identified only by positional coincidence with 
the pulsar's radio position (see Sect. ^). With about 
10 — 80 source counts the data of these pulsars do not 
allow a spectral analysis. Two other millisecond pulsars 
(PSR J0218-f4232 and J2124-3358) are detected by the 
ROSAT HRI only, which does not provide spectral in- 
formation. Nevertheless, an estimate for the pulsar's en- 
ergy flux and luminosity within 0.1—2.4 keV can be ob- 
tained from the observed count rate with only few ad- 
ditional assumptions: the pulsar's X-ray spectrum, its 
distance and the absorption column Nh- Detailed in- 
formation on the X-ray spectrum is only available for 
two (PSR J0437-4715, PSR 1821-24) of the nine de- 
tected millisecond pulsars. In both cases the best fitting 
spectral model implies a power-law nature for the emit- 
ted X-radiation. [Becker fc Triimpei] (1997) have shown 



recently that with the exception of the younger pulsars 
Crab, B1509— 58 and the Vela pulsar all pulse-phase av- 
eraged photon indices measured for rotation-powered pul- 
sars in the soft X-ray domain are consistent with a photon- 
index a w — 2. Calculating luminosities we therefore use 
the individual photon indices resulting from spectral fits 
and a canonical value of a = —2 where no spectral in- 
formation is available. The proper motion corrected pe- 
riod derivatives were used to compute E, B \ and the 



spin-down age r for PSR 1957 -1-20 (|Camilo et al. |1 994), 
PSR J0437-4715 ( |Bcll et al. |l995; ^andhu et al. [ 1997) 
and PSR J2124-3358 ( [Bailes ct a"i~| l997). 

The ROSAT PSPC and HRI temporal resolution of 
130 fin and 64 ^s, respectively, is sufficient to resolve X-ray 
pulses from millisecond pulsars. A low number of detected 
source photons does not necessarily preclude a search for 



(Becker et al. 1993a). For the reduction of the photon ar- 
rival times from the spacecraft coordinates and recorded 
times to the solar system barycenter and the barycen- 
tric dynamical time-scale (TDB), we performed the stan- 
dard procedures for ROSAT data (see Becker et al. 1993b) 



using the JPL DE200 Earth ephemeris and the pulsar 
ephemeris valid for the observational epoch. Correction of 
the photon arrival times for a pul sar's binary motion was 
performed using the method of Blandford fc Teukolsky 
(1976). 

Millisecond pulsars are stable clocks. Given the high 
precision of the pulsar's radio ephemeris, each X-ray 
photon arrival time can be directly related to the 
pulsar's rotation phase cj) using the relation (pi — 
fractional part of {fAti + ^ /Ai|) with i = 1 to TV, where 
N is the number of photons, / and / the pulsar frequency 
and its first time derivative at the reference epoch tref 
and At = ti — tref - The statistical significance for pulsa- 
tions was computed using the Z^-test with n = 1 to 10 
harmonics in combination with the H-Test to determine 
the optimal number of harmonics in the periodic signal 



(De Jager 1987; Buccheri fc De Jagei 1989). 



Various methods have been employed to measure the 
fraction of pulsed photons: approaches like fitting sinu- 
soids to the phase histogram using the Rayleigh-Power 
or estimating the minimum in a light curve to find the 
DC-level are often valid only for specific pulse shapes 
(e.g. real sinusoids) or completely neglect random fluctu- 
ations around the lowest level. An additional dependence 
on the number of bins used for the pulse phase histogram, 
makes the pulsed fraction often a rather method depen- 
dent estimate. To avoid these problems we have computed 
the pulsed fractions by using a bootstrap method. This ap- 
proach has been recently used by Swanepocl et al. ] (1996) 
in the analysis of gamma-ray pulsars. The advantage of 
this method is that the fraction of pulsed photons is cal- 
culated from the data without having to construct a his- 
togram or any other estimate of the light curve. 

Another concern is the number of phase bins used to 
construct pulse profile histograms for presentation pur- 
poses or for comparison with model light curves. Here, 
the bin width is limited on one side by the instruments 
temporal resolution but on the other side also by the sta- 
tistical significance of the signal. The latter is often not 
taken into account. In this work we have chosen the num- 
ber of phase bins to construct a pulse profile histogram by 
using the following approach: denoting the Fourier-power 
of the i-th harmonic by Ri and taking m as the optimal 
number of harmonics as deduced from the H-test, an exact 
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expression for the optimal number of phase bins is given 
by (De Jager 1998, in prep.) 



1/3 



M = 2.36 I ^i^Rf 



(1) 



This expression compromises between information lost 
due to binning (i.e. zero bin width to get all information), 
and the effect of fluctuations due to finite statistics per 
bin (i.e. bin width as large as possible to reduce the sta- 
tistical error per bin) . The total error (bias plus variance) 
is minimized at a bin width of 1/M. 

3. The faint detections 

3.1. The black-widow pulsar: PSR B1957+20 

The first millisecond pulsar reported in X-rays was the so 
called black-widow pulsar, a 1.6 ms pulsar at a distance 
of 1.53 kpc, which is in a close, approximately 9.16h or- 
bit with a low-mass white dwarf companion of ~ 0.025 
Mq (Fruchter et al. 1988). The source was observed in 
October 1991 with both the PSPC and the HRI detec- 
tor. We have reanalyzed the archival ROSAT data (Rev 
2) and found the source detected in the PSPC observa- 
tion with a total of 82 counts of which « 18 counts belong 
to the background. With an exposure time of 15 511s we 
deduce a PSPC count rate of (4.2 ± 0.6) x lO"^ ^ts/s in 
0.1—2.4 keV, approximately twice the number given by 
Kulkarni et al. (1992) for the range 0.5-2.0 keV. The pul- 
sar is marginally detected in the HRI data, in agreement 
with the results presented by Fruchter et al. [ (1992). An 
estimate of the pulsar's X-ray flux from the count rate 
results in fx — 3 X 10^^^ erg s^^cm^'^, implying a lumi- 
nosity of Lx 

X-ray efiiciency of L^/E « 0.8 x 10 
lar absorption we have adopted Nh 



8.5 X 103i((i/1.53 kpc)2 ergs-i and an 



For the interstel- 
4.5 X 10^1 cm-2 
as estimated by Kulkarni et al. | (1992) from the visual ex- 
tinc tion {Ay 1 — 2) of a star close to the pulsar position 



(see Pjorgovski fc Evans] 1988). 

To search for a modulation of the detected X-rays at 
the pulsar's rotation period, a photon arrival time analysis 
was applied to all photons selected from a 35 arcsec aper- 
ture centered on the pulsar's X-ray position. The selected 
counts represent more than 99.9% of the PSPC source 
photons with a background contribution of about 15%. 
Folding the corrected arrival times according to the pul- 
sar ephemeris valid for the ROSAT observational epoch 
did not reveal significant X-ray pulses. A not very restric- 
tive 1(T upper limit of 60% has been deduced for the X-ray 
pulsed fraction. 

Arons & Tavani93 (1993) have proposed that X-ray 
emission from PSR 1957-1-20 could arise from the pulsar 
itself or from the pulsar wind which is evaporating the 
companion star. The low number of counts recorded from 
this pulsar does not allow us to distinguish between the 
two mechanisms. 



3.2. PSR J0751+1807 

This binary pulsar has been recently discovered by 



Lundgren et al. (1995) in the error b ox of the uniden tified 
EGRET source CGRO J0749-M807 ( [Hartman et"al~| l992; 



Fichtel et al. | 1993) but the association between the pul- 



sar and the putative gamma-ray source could not be es- 
tablished (CGRO J0749-I-1807 is of marginal significance 
only and was finally removed fro m the 2EG and 3EG 
catalog of gam ma-ray sources ( Thompson et al. 1995; 
pJartman et al. 1998)). The 3.48 ms pulsar and its low- 
mass white dwarf companion are in a circular 6.3h orbit. 
The dispersion measure based pulsar distance is 2 kpc. 

ROSAT observed the pulsar in November 1993 
(13 992s, on-axis) and April 1994 r 7008s, 37 arcmin off- 
axis) with the PSPC ( Becker et al. | l996). Approximately 
50 source counts were recorded in the November 93 obser- 
vation leading to a PSPC count rate of (3.6 ± 0.6) x 10~^ 
cts/s. The degradation of the PSPC's point source sensi- 
tivity at 37 arcmin off-axis in combination with a reduced 
exposure time prevented a detection of the pulsar in the 
April 94 data. Estimating the pulsar's X-ray ffux and lu- 
minosity from the 1993 PSPC count rate yields fx = 8.3 x 
IQ-'^'^ erg s~icm^2 g^j^j = 4 x 10^^{d/2 kpc)^ erg s^^, 
implying an X-ray efficiency of L^/E m 5.3 x 10"'^. Con- 
verting the count rate to an energy flux is very sensitive 
to the assumed interstellar absorption. The pulsar's dis- 
persion measure implies a value of 9 x 10^°cm~^, mak- 
ing the usual but uncertain assumption that there are 
10 hydrogen atoms for each free electron along the line 
of sight. The model of Dickey & Lockman (1990) im- 



plies Nh ^ 5 X 10 cm for the total hydrogen col- 
umn throug h the Galaxy whereas from the HI survey of 
[Stark et al. | (1992) we deduce TV// ~ 4 x lO^" cm'^. The 
latter is adopted for the pulsar's absorption column. 

A timing analysis did not reveal signiflcant X-ray 
pulses at the expected radio period. A la upper limit of 
70% is deduced for the fraction of pulsed photons. 



3.3. PSR J 1012+ 5307 

PSR J1012+5307 is a 5.26 ms pulsar which is in a 14.5h 
binary orbit with a low-mass white dwarf companion. The 
pulsar was discovered recently during a survey for short- 
period pulsars using the Lovell radio telescope at Jodrell 



Bank (Nicastro et al. 1995). A detection of the companion 



star at optical wave-length was reported from an inspec- 
tion of the Palomar sky survey plates by the same authors. 
The distance of the pulsar as estimated from its dispersion 
measure is 0.52 kpc. 

A marginal X-ray detection of PSR J1012+5307 was 
reported recently by Halpcrn] (1996) who found the 
millisecond pulsar « 30' off-axis in a serendipitous PSPC 
observations of approximately 14.5 ksec duration. An ex- 
cess of « 80 ±20 counts above the background level within 
a circle of 2.5' and within 1' of the pulsar's radio posi- 
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tion has been interpreted as a possible detection of the 
pulsar (Halpern 1996). Converting the PSPC count rate 
to an energy flux we find fx ~ 4.9 x 10"^^ erg s^^cm^^ 
and = 1.6 x 10^°{d/0.52 kpc)^ erg s'^ for the X-ray 
luminosity within 0.1—2.4 keV. The latter implies an X- 
ray efficiency of L^/E « 0.4 x 10^"^. For the interstel- 
lar absorption we used Nh — 7 x 10^^ cm~^ as given by 
Snowden et al. (1994) for the HI hole in Ursa Major. No 
X-ray pulses have been detected from PSR J1012-I-5307. A 
Icr pulsed fraction upper limit is 75%. A 100 ksec ROSAT 
HRI observation to confirm the pulsars X-ray detection is 
scheduled for the second half of 1998. 



3.1 PSR J1024- 0719 and J1744-1134 

At present there are only nine galactic millisecond pulsars 
which are not in binary systems ( Camilo| 1996). PSR 
J1024-0719 and J1744-1134 ar e two of them . Both pul- 
sars were discovered recently by Bailes et al. (1997) dur- 
ing the Parkes 436 MHz survey of the southern sky. PSR 
J1024— 0719 has a rotation period of 5.16 ms and a period 
derivative of P = 1.84 x 10^^° s s^^, implying an upper 
limit to the pulsar age of > 4.4 x 10^ years and a rotational 
energy loss of > 5.3 x 10^'^ erg s~^. The spin parameters 
of PSR J1744-1134 are similar to those of J1024-0719: 
with a period of 4.07 ms and a period derivative of 
0.86 X 10^^° s s~^ the pulsars spin-down age is > 7 x 10^ 
years and the spin-down energy is > 1.42 x 10'^^ erg s^^. 
We note that for both pulsars the Shk lovskii contrib ution 
to P is not well constrained so far (see Bailes et al. 1997). 
The spin-down age and the rotational energy loss given 
above are preliminary for that reason. 

Both millisecond pulsars are at a relative close distance 
to us. For PSR J1024-0719 the Taylor & Cordes model 
implies d ^ 350 pc based on the pulsar's dispersion mea- 



sure dBailes et al. | 1997), whereas for PSR J1744-1134 the 
radio parallax has been measured recently and yielded a 
distance of about 260 ± 60 pc (Bailes 1997, priv. comm.). 

ROSAT HRI observations of PSR J1024-0719 and 
J1744— 1134 were performed in November/December and 
September 1997 for an effective exposure time of 80 635 
sec and 61 238 sec, respectively. The X-ray sources RX 
J1024.6-0719 (RA=10:24:38.60, DEC=-07:19:21.5) and 
RX J1744.4-1134 (RA=17:44:29.35, DEC=-11:34:52.2) 
were detected with a significance of ~ 4— Scr, only 2—3 arc- 



sec distant from the pulsars' positions (Becker & Triimper 
1998; [Becker et ar| l998, see Fig.0). The positional off- 
set between RX J1024.6-0719, RX J1744.4-1134 and 
the pulsars' radio position is well within the uncertainty 
of the satellite pointing accuracy. The probability that 
RX J1024.6-0719 and RX J1744.4-1134 are chance su- 
perpositions of unrelated background object is ^ 10~^ 
and ^ 3 x 10^^, respectively, estimated from the density 
of X-ray sources in the HRI field of view. We therefore 
accept RX J1024.6-0719 and RX J1744.4-1134 as the 




Fig. 1. Sub-image of the PSR J1024-0719 field as ob- 
served with the ROSAT HRI in November/December 
1997. The image has a spatial binning of 5 arcsec. The pul- 
sar's X-ray counterpart RX J1024.6— 0719 is indicated by 
an arrow. Several unidentified X-ray sources are detected 
in the pulsar's neighborhood, only few arcmin distant from 
the pulsar 

likely X-ray counterparts of the millisecond pulsars PSR 
J1024-0719 and J1744-1134. 

Converting the HRI source count rate of (3.3 ± 
0.9) X lO-'' and (2.9 ± 0.9) x 10-^ as deduced for 
J1024-0719 and J1744-1134, to an energy flux we 
find - 2 X lO-i"* ergs-^cm-2 within 0.1-2.4 
keV for both sources. The latter implies an X-ray lu- 
minosity of ij; — 3 X 10^^ (d/0.35 kpc)^ erg s~^ and 

= 2 X 10^9 (d/0.26kpc)2 erg s'^ for J1024-0719 and 
J1744— 1134. For the interstellar absorption we used Nh = 
2 X 10^° cm~^ and 10^° cm~^ respectively, as deduced from 
the pulsars' dispersion measure. The number of recorded 
counts was not sufficient to perform a search for pulsa- 
tions. 



4. The nearby and bright pulsar PSR J0437-4715 

The first millisecond pulsar for which a more detailed 
knowledge of its X-ray emission properties has be- 
come available is the 5.75 ms pulsar PSR J0437— 4715. 
The source, which was discovered in the Parkes south- 
ern sky survey by Johnston et al. 
is in a close 5.74-day circular 
0.25 Mq white-dwarf companion ( Johnston et al. [ 1993; 

Danziger et al~[L993) 



(1993) and which 
orbit around a « 



Becker et al. 1993a; Bailyn 1993 



is the nearest and brightest millisecond pulsar known. 
Recent measurements of the annual parallax have 
adjusted the pulsar's distance to be 178 ± 26 pc 
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( ^andhu et al. 1997), placing it somewhat further away 
than originally indicated by the dispersion measure and 



the electron density model of Taylor & Cordes (1993). The 
parallax distance further constrains the pulsar's proper 
motion and revises its intrinsic period derivative down 
to Pint — (0.8 ± 0.7) X 10^^", approximately a factor 
of 7 smaller than the observed one Pobs = 5.73 x 10"^" 
(Bell ct al. 1997). Based on Pint the pulsar's spin-down 



energy and polar magnetic field is found to be log£' = 
33. 22;° J erg/s and \ogB±_ = 8.34^°;;^!^ G, respectively. 
The pulsar's characteristic age is P/2Pint < 6 x 10^ yr. 

X-ray emission from the pulsar was first detected in 
the ROSAT all-sky survey by [Becker et al. | (1993a). Based 
on a serendipitous sa 6 ksec pointed observation with the 
ROSAT PSPC [Becker fc Triimpcij (1993) found that the 
X-ray emission is pulsed and the pulsar's soft X-ray spec- 
trum is best described by a single power-law. An energy 
dependent pulsed fraction showing a peak at ^ 0.9 keV 
was found by fitting a sine wave to the X-ray pulse profiles 
obtained for different energy ranges. Because it appeared 
difficult to explain this apparent energy dependence in 



terms of single-component emission Becker & Triimper 



(1993) proposed that the radiation from PSR 0437-4715 
is a mixture of spin-modulated thermal emission from 
heated polar caps combined with unpulsed emission from 
a pulsar wind or plerion. A reanalysis of the serendipitous 
ROSAT data together with data tak en with the EUVE 
satellite at 65 — 120 A (Halpern et al. 1996) confirmed the 



results of Becker & Triimper (1993) but has weakened the 



case of an energy dependence of the pulsed fraction by us- 
ing improved pulsar ephemeris and a different method to 
estimate the fraction of pulsed photons. 

In the following we present the results from a re- 
analysis of the archival ROSAT data of PSR 0437-4715 
including unpublished data from subsequent observations 
with both the ROSAT PSPC and HRI detectors. 

4.1. Observations 

Since its discovery in 1993, PSR J0437-4715 has been the 
target of several pointings with both the PSPC and the 
HRL A total of - 20 ksec PSPC observations and - 60 
ksec HRI observations have been performed with different 
sensitivities. The coverage of the pulsar's binary orbit by 
ROSAT observations is illustrated in Fig. |^. 

Pointed observations with the ROSAT PSPC were per- 
formed on 1992, September 20—21 for an effective expo- 
sure time of 5846 s and on 1994 July 2-4 for 9885 s and 
4459 s, respectively. The latter observation was done with 
the boron filter. As a consequence, the pulsar's PSPC 
count rate is reduced to 0.042 ± 0.003 cts/s, compared 
with a count rate of 0.204 ± 0.005 cts/s for the filter-off 
observation. In addition, ROSAT observations with the 
HRI were performed between 1994 July 18th and August 
24th and on 1997 September 22—26 for an effective expo- 
sure time of 35 209 and 22 928 seconds, respectively. In the 
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Fig. 2. Coverage of the pulsar's binary orbit by X-ray 
observations using ROSAT. Circles indicate observations 
with the PSPC in the focus of the telescope whereas tri- 
angles represent HRI observations. 



later observation the satellite wobble-mode was disabled. 
The pulsar's HRI source count rate is 0.042 ± 0.002 cts/s. 

4-.2. Spatial emission properties of PSR J04^37^4715 

PSR J0437— 4715 is one of few millisecond pulsars for 
which an Hn bow-shock nebula has been detected 



(Bell et al. 1993; 199 5). Others are the black- widow pul- 
sar PSR 1957+20 ( [Kulkarni fc Hestei) 1989) and PSR 
J2225+ 6535 which is s urroi mded by the s o called Guitar 
Nebula ( Cordes et al. 1993; Romani et al. 1997). It is gen- 
erally believed that the pulsar wind, a mixture of charged 
particles and electro magnetic radiati on, shocks the inter- 
stellar medium (see Arons fc Tavan: 1993) giving rise to 
a bow-shock. Such a shock region can produce both ther- 
mal emission and synchrotron radiation. Only the local 
conditions, such as the magnetic field and particle den- 
sity, determines if X-rays are produced and which emission 
process dominates. 

The optical bow-shock nebula of PSR J0437-4715 has 
a standoff distance from the pulsar position of about 12 
arcsec (see Fig.||). Diffuse optical emission of the nebula is 
visible up to about 1 arcmin from the pulsar position along 
the side wings of the bow-shock hypcrboloid. The pulsar's 
proper motion direction is 122° ± 2° ( [Bch ct al. [ l995), 
measured anticlockwise from north. 

With a spatial resolution of ~ 25 arcsec (FWHM) PSR 
J0437— 4715 appears to be point like in the available PSPC 
data. A source extension at energies below ~ 0.2 keV is 
found to be caused by the well known ghost imaging effect. 
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Fig. 3. Sub-image of the pulsar field as observed with 
the ROSAT HRI in September 1997 (satellite wobble dis- 
abled). The image has a spatial binning of 2.5 arcsec. The 
4.3 arcmin neighboring source RX J0438. 1-4710, identi- 
fied to be a Seyfert 1 Galaxy at z = 0.051 ( |Grupc 1996) 
is visible to the upper left. The inset in the lower right 
corner shows the pulsar's optical bow-shock nebula as ob- 
served by A. Fruchter in Ha- The arrow indicates the pul- 
sars proper motion direction, which has an orientation of 
122° measured anticlockwise from north. 



an instrumental artefact which mimics an extended source 
for ultra-soft photons. Compared with the ROSAT PSPC 
the spatial resolution of the HRI detector is higher by 
about a factor of 5. The HRI therefore is more suitable 
to search for an extended emission component on small 
angular scales. Inspecting the 1994 HRI data we found a 
small elongation of the pulsar's X-ray counterpart towards 
a direction of 105± 5° measured anticlockwise from north. 
Since this orientation is consistent with the satellites wob- 
ble direction 102°) we conclude that this elongation is 
an artefact of uncorrected residual wobble rather than a 
feature associated with the pulsar's Ha nebula. To pre- 
vent such a source extension due to incomplete attitude 
correction we observed the pulsar in 1997 with the satellite 
wobble disabled. Figure || shows the HRI image of the pul- 
sar field based on the 1997 data. As can be seen from this 
image there is no indication of an ellipsoidal elongation in 
the orientation of the optical bow-shock nebula as found in 
the 1994 HRI data. Figure ^ shows that the source extent 
of the pulsar's X-ray counterpart is in agree ment w ith the 
ROSAT HRI point-spread- function (see also Boese 1998). 
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Fig. 4. a. The radial surface brightness distribution 
within an 80 arcsec circle centered on the pulsar PSR 
J0437-4715 and based on the 1997 HRI data. The his- 
togram incorporates a point source and an additional 
broadening of the PSF to account for an uncorrected 2.6 
arcsec attitude error, smearing out the model point spread 
function, b. Residuals of the PSF fit in units of a. The fit 
has a reduced of 0.87. 



To search for a possible X-ray emission component 
from the side wings of the bow-shock nebula we used 
the merged 1994 and 1997 HRI data. 785 counts were se- 
lected from within three quarter of a ring around PSR 
J0437— 4715, having an inner radius of 12" and and outer 
radius of 60" in the orientation of the optical image of 
the pulsar bow-shock. Figure || shows the selected sky 
region scaled up by a factor 6 compared with Fig. |[ 
Its area is approximately 2.6 arcmin^, implying a surface 
brightness of ^ 302 cts arcmin"^. Shifting the mask to 
source free regions in the neighborhood of PSR J0437- 
4715 we find that on average 642 of the 785 counts in the 
mask belong to the background. The excess of 143 counts 
above the background level is about 6% of the ~ 2480 
pulsar counts falling within the 12 arcsec inner radius of 
the mask and thus is well explained by mirror scattering 



(David et al. 1997). For the 3a count rate upper limit of a 
diffuse nebula component we find 3 x -\/785/58 138 cts/s = 
1.4x 10^'^ cts/s (see Becker 1995, p65). Assuming a power- 
law spectrum with a photon-index a = —2 and a col- 



umn density of Nh = 0.8 x 10^°cm~^ (see Sect. 4.3) 



the count rate upper limit corresponds a 3a flux up- 
per limit of < 5 X 10"^'' erg s~^ cm~^ and a luminosity 
of < 2 X 10^9 erg s-i (d/180pc)2 in 0.1-2.4 keV. Com- 
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Fig. 5. Sky region around PSR J0437-4715 selected to 
quantify a possible X-ray emission component from the 
side wings of the pulsar's bow-shock nebula. The mask 
has an inner radius of 12" and and outer radius of 60" 
in the orientation of the optical image of the pulsar bow- 
shock. The pixel size is 5x5 arcsec. Based on the merged 
1994 and 1997 HRI data we find in total 785 counts within 
the covered area. 

pared with the pulsar's spin-down energy this is less than 
0.2 X 10~^ E, similar to the 3cr conversion upper limit 
found for a possible bow-shock contribution to the X-ray 



emission from 1957-1-20 (Fruchter et al. 1992) 



4-.3. Spectral emission properties 

Becker & Triimper (1993) have shown that the pulsar's 
soft X-ray spectrum in 0.1—2.4 keV is best described 
by a power-law dN/dE cx i?". Figure || shows the 
contour plot of this model in the Nh — a plane. The 
best fitting parameters for this model, deduced from the 
merged 1992 and 1994 PSPC data, are a column den- 
sity Nh = (8 ± 2.5) X lO^^cm"^ and a photon-index 
of a = —2.35 ± 0.15. This is in agreement with the re- 
sults foimd by Halpern et al. (1996) based on a combined 
EUVE/ROSAT spectral analysis using the Sep. 1992 data 
only The 1994 PSPC data taken with the Boron filteiQ 
were used for a cross check of the fitted power-law spec- 
tral parameters by computing the count rate ratio for an 
observation with and without a Boron filter and compar- 
ing this with the measured value. The result was found to 
be in agreement with the filter-off spectral fittings. 

The pulsar's energy flux deduced from the best fitting 
power-law model is fx — 1.9 ± 0.2 x 10^^^ erg s^^cm^^ 
in 0.1—2.4 keV, implying an X-ray luminosity of Lx — 
7.3 ±0.8 X 10^° erg s~^ (d/180pc)^ and a conversion factor 
of L^/i;^4x 10-3. 

It was already shown by [Becker fc Triimpei (1993) 
that a single black-body spectrum does not fit the ROSAT 
data and leaves a residual hard excess above ~ 0.5 keV 



^ The advantage of the Boron filter is to improve the PSPC's 
spectral resolution in the carbon band (0.1 — 0.288 keV) by a 
further subdivi sion of this energy range into two bands. See 
Stephan et al. L991 for informations on the filter transmission. 



PSR J0437-4715 




Fig. 6. Contour plot of as a function of the Galactic 
column density Nh x 10^^ cm"^ and the photon-index a. 
The confidence levels range from la to bigma. 



see also Becker 1995, pl25), whereas a double black-body 



model yiel ds acceptable spectral fits. Based on ASCA ob- 
servations, Kawai et al. (1998 - 



hereafter KTS98) have re- 
ported recently that the pulsar spectrum beyond 0.7 keV 
fits better with a black-body rather than with a compos- 
ite (power-law plus black-body) or simple power-law spec- 
trum. The authors deduced a temperature of T ^ (3 ± 
0.5) X 10^ K, a black body radius of Rbb - 54 m (d/180 pc) 
and a bolometric flux of fbb ~ 5.1 x lO-^'^erg s"^ cm"^. 
Figure ^ shows the ROSAT PSPC data fitted with a dou- 
ble black-body spectrum in which the parameters of the 
second (harder) spectral component have been fixed ac- 
cording to the ASCA results. The parameters of the soft 
black-body component have been fitted and imply a tem- 
perature of 1.2 X 10^ K, i?;,6 ~ 510 m (d/180pc) and a bolo- 
metric flux of /f,;, — 8.6 X 10"^^ erg s"^ cm^^. The column 
density in this model is fitted to be zero. Due to the small 
radii both black-body components are attributed to the 
pulsar's polar cap. The X-ray flux detected with ASCA 
above 0.7 keV is not pulsed. The authors give a pulsed 
fraction upper limit of 70%. 

An alternative to the double black-body model, also 
suggested by KTS98, is that the simple power-law com- 
ponent seen by ROSAT and EUVE has a spectral break 
near the upper boundary of the ROSAT band. In fact, 
the authors found that both a broken power-law model 
or a power-law model with exponential cut-off can de- 
scribe the X-ray data from PSR J0437— 4715 very well. 
The ASCA data imply a spectral break (cut-off) near 
'--^ 1.4 keV with a slope of about —4 for the second power- 
law component. The broken power-law model fitted to the 
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Fig. 7. a. The ROSAT PSPC data of PSR J0437-4715 fit- 
ted with a two component black body model. The harder 
component (dashed line) represents the black-body spec- 
trum as deduced from spectral fits to ASCA data (see 
text). The dotted line represents the soft component as 
obtained from a spectral fit to the ROSAT data. b. The 
PSPC data fitted with a broken power-law. The break 
energy in this model is 1.4 keV. The slopes of the first 
and second power-law components are —2.35 and —4.5, 
respectively. 

ROSAT data is shown in Fig. 0b. The slopes of the first 
and second power-law components are —2.35 ± 0.15 and 
—4.5 ± 1.65, respectively. The column density is found to 
be Nh = 0.8 X 10^0 cm-2. 

A number of different spectral models which assume 
thermal emission from polar hot spots taking into ac- 
count the influence of a possible neutron star H-, He- 
or Iron-atmospheres have been tes ted recently for PSR 
J0437— 4715 by several authors (see Rajagopal fc Romani 



Data 


Range 


Source 


Bg 


Pulsed 


PF 


set 


[keV] 


cts 


cts 


cts 


[%] 




0.1 - 2.4 


1200 


24 


326 ± 48 


28 ± 4 


PSPC 


0.1 - 0.6 


896 


19 


232 ± 38 


27 ± 4 


1992 


0.6 - 1.1 


211 


3 


67± 19 


32 ± 9 




1.1 - 2.4 


102 


2 


16± 14 


18 ± 14 




0.1 - 2.4 


2001 


80 


588 ± 74 


29 ± 4 


PSPC 


0.1 - 0.6 


1594 


67 


410 ± 59 


27 ± 4 


1994 


0.6 - 1.1 


363 


10 


107 ± 33 


30 ± 9 




1.1 - 2.4 


158 


8 


27± 20 


16 ± 14 


Boron 


0.1 - 2.4 


145 


4 


37± 13 


26 ± 9 


HRI-94* 


0.1 - 2.4 


1055 


9 


300 ± 50 


29 ± 5 


HRI-97 


0.1 - 2.4 


985 


12 


256 ± 40 


26 ± 4 



Errors represent the la confidence range. 
'Based on data taken on 21. and 22. Aug. 1994 



1996; [Pavlov et ar| l996; [Zavlin fc Pavlov) 1998). All these 
models yield spectral fits which arc in a somewhat better 
agreement with the data than a simple black-body spec- 
trum but still leave some unmodeled hard excess above 1 
keV. The polar cap model of Zavlin fc Pavlov| (1998) pre- 
dicts a pulsed fraction of more than 50% for the emission 
above 1.1 keV, in contradiction to what is observed with 



ROSAT (see Sect. 14 and Table 



4-. 4- Temporal emission properties 

Because of the long time gap between the two PSPC ob- 
servations and the difference in detector type between the 
HRI and PSPC we have analyzed the available ROSAT 
data from PSR J0437— 4715 independently of each other. 
The results can therefore be regarded as independent mea- 
surements. In the following we give a detailed description 
of the data analysis for the different data sets. 

Quantities such as the number of counts, the back- 
ground contribution and the number of pulsed counts as 
well as the fraction of pulsed photons for different energy 
ranges are summarized in Table |l|. 

4.4.1. The 1992 PSPC data 

A timing analysis based on this data but using a pre- 
liminary pulsar ephemeris led to the discovery of the X- 
ray pulses from PSR J0437— 4715 by Becker & Triimpei 
(1993). The results of a reanalysis were presented recently 
by Halpern et al. (1996) using an improved ephemeris. 



However, to perform a timing analysis consistently for all 
five data sets we have analyzed the PSPC-92 data again 
using a pulsar ephemeris valid at an epoch close to the 
ROSAT observational epoch. 



W. Becker & J. Triimper: The X-ray emission properties of millisecond pulsars 



9 



0.1-0.6 keV 



PSPC-92 



140 



120 - 



J2 100 - 




1.1-2.4 keV 



35 



25 - 



20 

15 ■ 

40 - 

30 



c) 



0.1-2.4 keV 
Boron 



10 -'- 



d) 



0.0 0.2 0.4 



0.6 0.8 1 

Pulse 



1.2 

Phase 



1.6 1.8 2.0 



200 



_ 0.1-2.4 keV 



150 - 



100 



50 



3 150 



100 



- f) 



200 
150 - 
100 

1.4 

1.2 
1.0 

0.8 
0.6 
0.4 
0.2 

0.0 




e) 



HRI- 



H 1 h 



0.1-2.4 keV 





0.1-2.4 keV 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1 
Pulse Phase 



.8 2.0 



Fig. 8. Integrated pulse profiles of PSR J0437-4715 for different energy ranges. Two phase cycles are shown for clarity 
and la errors are indicated. Panels a-c depict the profiles based on the 1992-PSPC data. Panels d-g show the pulse 
profiles within the full ROSAT energy range taken with the HRI in 1994 (e) and the PSPC in 1992 and 1994 (d,f,g), 
respectively. The DC-level and its la error is indicated by solid and dashed lines, respectively. No phase shift is observed 
for the pulses within the ROSAT energy band. The X-ray pulses are found to be consistent with a constant pulsed 
fraction and a pulse shape which appears similar to the radio one shown in panel h (see text and Table ^ for details) . 
The relative phase between the X-ray and radio profiles in this figure has been fixed arbitrarily (see Sect.|X|). 
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For the photon arrival time analysis of the 1992-PSPC 
data we have selected 1200 photons from an 70 arcsec 
aperture centered on the pulsar's X-ray position. The cut- 
radius was chosen according to the PSPC's point-spread 
function and contains more than 99.9% of the pulsar 
counts. 70 arcsec is large enough to select most of the 
ultra-soft photons subject to the electronic ghost-imaging 



effect (Bricl ct al. 1996) which is barely present in this 
data set. An extraction radius of 100 arcsec as used 



by Halpern et al. (1996) appears too large and selects an 
unnecessary large fraction of background photons. 

The folded pulse profiles for different energy ranges are 
shown in Fig. ||a— c. The H-Test yields the highest proba- 
bility for 2 harmonics, indicating a small deviation from a 
sinusoidal shape. The width of the X-ray pulse is approx- 
imately 0.8 rotations (« 290°), practically the same as 
observed in the radio domain (see Fig.^!) . The fraction of 
pulsed photons, measured for the different energy ranges 
by using the bootstrap approach is given in Table ^ There 
is no significant pulsed X-ray emission above 1.1 keV in 
the 1992-PSPC data (see Fig.|c). Only 102 photons have 
been detected in this energy range, for which the boot- 
strap method indicates a pulsed fraction of 18 zt 14 %. In 
contrast to our earlier finding s ([Becker fc Triimpei| 1993) 



the bootstrap method indicates that the pulsed fraction is 
independent of the photon energy in the ROSAT band. 

4.4.2. The 1994 PSPC data 

In contrast to the 1992-PSPC data, we find that the spa- 
tial extent of the pulsar's X-ray counterpart below ~ 0.25 
keV is strongly dominated by the electronic ghost-imaging 
effect. In order not to exclude these ultra-soft photons 
from the analysis we have selected the pulsar counts below 
0.25 keV from an extended region chosen according to the 



expected pattern of photon imaging (Briel et al. 1996). 
The PSPC count rate from this area yields a rate of 
0.205 ± 0.005 cts/s, in agreeme nt with the count ra te ob- 
served in the 1992-PSPC data ([Becker fc Triimper| 1993). 
In total, we have selected 2101 counts for the timing anal- 
ysis, more than 99% of the total pulsar contribution. Fold- 
ing the corrected and barycentered arrival times with the 
pulsar's rotation frequency resulted in the pulse profile 
shown in Fig.||g. The pulse profiles obtained for the en- 
ergy ranges 0.1-0.6 keV, 0.6-1.1 keV and 1.1-2.4 keV 
are not significantly different from the profiles shown in 
panel ^a-c. The pulse profile obtained from the observa- 
tion using the boron filter is shown in Fig. [^d. 

4.4.3. The 1994 and 1997 HRI data 



As already pointed out in Sect. 4.1, the pulsar's 1994 HRI 
observation of ~ 35 000 s is distributed over a period of 
38 days. Approximately 9800 s of the exposure time was 
taken between July 18th and August 13th 1994 whereas 
about 25 300 s of the exposure was taken within only two 



days on August 21-22. Selecting the photons from an aper- 
ture of 10 arcsec radius centered on the pulsar yielded 1484 
counts available for the timing analysis. These are approx- 
imately 90% of all pulsar counts and a background con- 
tribution of only « 1%. Folding the arrival times with the 
pulsar's rotational frequency (using the same ephemeris 
as already used for the timing analysis of the 1994-PSPC 
data) we find from the H-Test that the significance for pul- 
sations is maximized for 3 harmonics at « 91, indicat- 
ing a somewhat stronger deviation from a sinusoidal pulse 
shape than previously found in the PSPC data. However, 
the fraction of pulsed photons is found to have a rather 
low value of only 20 ±4%. For a consistency check we have 
therefore performed the same analysis again but at this 
time restricted to the data taken on 1994 August 21. and 
22. only. The result is interesting in the sense that we find 
both the significance of the pulsed signal and the fraction 
of pulsed photons enhanced to Z| 120 and 29 ± 5%, re- 
spectively, although the number of source counts available 
for the analysis is reduced to 1055. The shape of the pulse 
profile found in both data sets is not significantly differ- 
ent and shows only a change in the DC level. The profile 
based on the 21—22 August 94 data is shown in Fig.^. 

The differences in significance and pulsed fraction are 
probably due to the limited stability of the ROSAT on- 
board clock having a nominal drift of w 8 x 10~* per 
day. At the pulsar's observational epoch, clock calibration 
measurements were only performed once per week. Inves- 
tigating the residuals of a linear fit from the space craft 
clock calibration against UTC we find a change in the drift 
rate approximately in the middle of the first half of the 
1994 HRI observation whereas the drift rate stays nearly 
constant for the remaining time of the observation. From 
this we conclude that folding all data from the full 38 
day observation span results in a smearing of the pulsed 
emission leading to a reduction of the pulsed fraction. We 
conclude that the pulsed fraction of 29 ± 5% for the 1994 
HRI data is in agreement with the results of the two PSPC 
observations. 

For the HRI data taken in September 1997 we applied 
the same analysis as for the 1994 data. The results are 
in agreement with findings deduced from the 1994 HRI 
observation (see Table |l|). 

4.4.4. Relative phase between the X-ray and radio pulse 

In view of the short pulse period of millisecond pulsars 
and the uncertainty in the ROSAT space craft clock cal- 
ibration against Coordinated Universal Time (UTC), a 
comparison of the relative phase between the X-ray and 
radio pulses for PSR J0437— 4715 appears difficult. Cali- 
bration measurements of the satellite's space craft clock 
against UTC are available with a frequency of usually once 
per week for the pulsar observational epochs in 1992 and 
1994. A daily clock calibration has been performed only 
since November 1996. 
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In order to measure the X-ray pulse arrival relative to 
the radio pulse we made local fits to the ROSAT clock 
calibration points against UTC using three and four order 
polynomial functions. With this approach we find resid- 
uals which are in the order of 1 — 2 ms. Although this 
surely represents the limits of the clock accuracy against 
UTC it still corresponds to a relative error in the X-ray 
pulse arrival time of approximately 20 — 30%(lcr). In view 
of these uncertainties the relative phase between the radio 
and X-ray pulse of PSR J0437-4715 cannot be constrained 
with the available ROSAT data. More accurate timing will 
be available on AXAF and XMM, which will provide the 
desired information. 



5. The globular cluster pulsar B1821-24 in M28 

PSR B1821— 24 is an isolated 3.1 ms pulsar in the glob- 
ular cluster M28. Among the few tens of known millisec- 
ond pulsars it has the highest spin-down energy. X-ray 
emission from the source was detected with the ROSAT 



PSPC in March 1991 (|Danncr ct al. | l994) and with the 
ROSAT HRI and the ASCA CIS & SIS detectors in March 



1995 dPanner et al. | 1997; |Saito et al. | 1997). The results of 
these observations are remarkable in the sense that they 
detect pulsed X-rays in the wide energy band from 0.1 — 10 
keV. 

Figure ^ shows the X-ray and radio pulse profile as 
observed with the ROSAT HRI (top profile), the ASCA 
CIS-detector (middle profile) and the National Radio As- 
tronomy Observatory in Creen Bank (bottom profile) . The 
X-ray light curves are double peaked with a very narrow 
profile for the pulse components. The significance of the 
X-ray pulses is about Aa in the HRI data. 

The high resolution puls e profile observed a t 800 MHz 
was published recently by Backer & Sallmen (1997). At 
this frequency, the radio profile shows three pulse compo- 
nents, of which peak #2 is found to show intensity vari- 



ations on time scales of hours to days. Backer & Sallmen 



(1997) found that the spectral-index of the third and 
broad radio pulse component is similar to the one ob- 
served for the narrow and leading pulse component #1. 
The spectral-indices measured between 470 - 1700 MHz 
are —1.4 ± 0.2 and —3.1 ± 0.2 for the first and second 



pulse component (Foster et al. 1991). 

Unfortunately neither the ROSAT nor the ASCA clock 
calibration is sufficient to phase relate the X-ray pulses of 
a 3.1 ms pulsar with the radio pulse. However, recent ob- 
servations of PSR 1821—24 with the Rossi X-ray Timing 



Explorer (Rots et al. 1997) have shown that the primary 
X-ray pulse component appears to be nearly phase aligned 
with the radio pulse component dominating at 800 MHz. 
Phase alignment between radio and X-ray pulses is only 
established (with sufficient accuracy) for the Crab pulsar. 
Unless the phase alignment is chance coincidence it sug- 
gests a common emission site for both the main radio and 
the dominating X-ray pulse. 
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Fig. 9. Integrated pulse profiles of the globular cluster 
pulsar P SR 1821-24 a s observed with the ROSAT HRI 
(top) by [Danncr et al. | (1997), the ASCA CIS detector by 
[Saito et al. |(1997) and with the NRAO at 800 MHz (bot- 
tom) by Backer fc Sallmei^ (1997). Two phase cycles are 
shown for clarity. The X-ray pulse profiles are character- 
ized by a double peak structure with a phase separation 
of ~ 0.6 between the two peaks. The radio profile at 800 
MHz depicts three pulse components. At this frequency 
the dominating radio pulse is nearly phase aligned with 
the primary X-ray pulse. 



The spectral analysis based on the ASCA data shows 
that the pulsars' X-ray spectrum between 0.7 — 10 keV 
is well described by a power-law with a phase averaged 
photon-index of about a ~ —2 ( ^aito et al. 1997). The 
first and the second pulse show photon- indices of — 1.5±0.3 
and —2 ± 0.3, respectively. The spectrum of the globular 
cluster emission, which is assumed to be the interpulse 
emission, is not well constrained. The interpulse emission 
appears to be softer than the pulsar emission, but the data 
do not allow us to discriminate between several possible 
spectral models: a power-law with a photon-index of —2, 
an optically thin thermal model with kT=1.4 keV and 
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a kT=6.2 keV bremsstrahlung model yield all acceptable 
fits. 

Imaging M28 with the ROSAT HRI has shown two 
separate sources: the point source RX J1824.5-2452P con- 
sistent with the millisecond pulsar position and a brighter 
extended source RX J1824.5-2452E, whose nature is not 



yet clear (Danner et al. 1997). The presence of the two ob- 
jects, unresolved by ASCA, further confuses the spectral 
interpretation. The power-law spectrum of the interpulse 
data (i.e. the spectrum of the extended source RX J1824.5- 
2452E) favors the model in which RX J1824.5-2452E is a 
pulsar powered synchrotron nebula, similar to the Crab. 
An alternative and may be more likely interpretation, 
however, is that RX J1824.5-2452E is made of a number 
of point sources (e.g. accreting binaries containing white 
dwarfs or neutron stars) which could not been spatially 
resolved by the HRI. An interpretation in terms of low ac- 
cretion LMXBs is also supported from a long term X-ray 
luminosity study of M28. Our reanalysis of the archival 
PSPC and HRI data, taken in March 1991 and 1995, sug- 
gest an X-ray flux variability of M28 on time scales of 
years. The total ROSAT HRI energy flux taken in March 
1995 from M28 is about a factor of 3 higher than the total 
energy flux deduced from the March 1991 PSPC obser- 
vation. This behavior is in line with the results recently 



published by Gotthelfc & Kulkarni (1997) who discovered 



an X-ray burst from M28 in the 1995 ASCA data. 

Estimating the pulsar's energy flux and luminos- 
ity from the HRI count rate (3.7 ± 0.4 x IQ-^cts/s) 
yields for the power-law spectrum fx = (5.2 ± 
0.6) X 10^13 erg s-^cm-2 and = 1.6 ± 0.2 x 

10'^'^(d/5.1 kpc)^ erg within 0.1—2.4 keV, respec- 
tively. The latter implies an X-ray efficiency of L^/E w 
0.8 X lO"'^. For the interstellar absorption we used Nh = 
2.8 X 10^^ cm~^ which is consistent with the ASCA mea- 
surements dSaito ct al. 1997) and the column density de- 
duced from the dispersion measure. 



6. PSR J0218+4232 



The 2.3 ms pu lsar PSR J0218 
Navarro et al. J(1995) 



4232 was discovered by 
The pulsar is in a two day binary 



orbit with a low-mass white dwarf companion and shows 
significant unpulsed radio emission throughout the pulse 
period. The latter has been taken by Navarro et al. (1995) 



as an indication that the magnetic dipole is almost aligned 
with the rotation axis. The Taylor & Cordes model indi- 
cates a dispersion measure based distance of 5.7 kpc. 
X-rays from the pulsar were detected on the basis of 



ROSAT HRI observations by [Verbunt et al. | (1996). Our 
re-analysis of the archival ROSAT data confirms their 
findings. Approximately 47 source counts were recorded 
from PSR J0218-h4232 during the 22 035s HRI observation 
in August 1995, implying a count rate of (2.1±0.4) x 10"^ 
cts/s. This HRI observation was performed 12 arcmin 
off-axis, resulting in a reduced spatial resolution and point 
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Fig. 10. The X-ray and radio pulse profiles of PSR 
0218+4232 as observed with the ROSAT HRI (top) by 
pCuiper et al. (1998) and with the Effelsberg radio tele- 
scope at 1409 MHz (bottom) by Lange et al. (1998, in 
prep.). Two phase cycles are shown for clarity. The relative 
phase between the X-ray and the radio pulse is unknown. 
In this plot it has been arbitrarily fixed. 

source sensitivity. The source is also marginally detected 
in a serendipitous ~ 25 ksec PSPC observation near the 
edge of the detector's field of view. 

X-ray pulses from PSR J0218-h4232 have been de- 
tected only recently using data from a subsequent 
ROS AT HRI observ ation which was performed in July 
1997 ( Kuipcr et al. 1998). About 182 source counts were 
recorded during an exposure time of 98096 s which re- 
vealed the existence of a pulsed signal with a statistical 
significance of ~ 5cr and a pulsed fraction of 37± 13%. The 
integrated pulse profiles observed between 0.1-2.4 keV and 
in the radio channel at 1409 MHz are shown in Fig. ^ 
As for PSR 2124-3358 (see Sect. 0) there is some indica- 
tion for the existence of two pulse components in the X- 



ray lightcurve. Kuiper et al. (1998) mention that the 1998 
ROSAT HRI data reveal some marginal evidence for the 
existence of a compact synchrotron nebula surrounding 
the pulsar. The authors found in their analysis a source 
extent of about 14 arcsec. However, the H RI point spread 
function (PSF) assumed by Kuiper et al. (1998) may have 
been underestimated by neglecting systematic uncertain- 
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ties in the HRI-PSF due to incomplete attitude correction 
(see Boesc| 1998 and discussion therein). 

No information on the pulsars X-ray spectrum is avail- 
able so far. The HRI provides no spectral information 
and the number of counts recorded in the serendipitous 
PSPC observation foes not allow spectral modeling. An 
estimate of the energy flux and luminosity obtained from 
the HRI count rate yields fx = 1.5 x 10""'^^ erg s~^cm" 
and Lx = 5.7 x 10^2^^/5 7 ^.^^^2 grg s^^ in 0.1-2.4 
keV, implying an X-ray efficiency of L^/E « 2.3 x 10"'^. 
The pulsed energy flux is found to be = (3.9 ± 

1.4) x 10^^* erg s~^cm" " 



stellar absorption was adopted from Verbunt et al 
who 



dKuiper et al. |1998). The inter- 

(1996) 



used the maximum color excess E{B — V) = 
0.09 dBurstein fc Hcilcs| 1982) in the direction of PSR 
J0218+4232;^v ^ ' S .lEjB-V) and Nh = 1.79x10^1 Ay 
(iPredehl fc Schmitt] 1995) which results in Nh « 0.5 x 



10^1 cm-^. 



Although the second EGRET catalog lists a source 
whose 8 0% confidence co ntour contains the m illisecond 
pulsar ( Thompson et al. 1995; Verbunt ct al. 1996), the 
third EGRET catalog ( [Hartman et al. | 1998) which is 
based on more viewing periods and better data statistics 
than the 2EG catalog invalidates this putative identifica- 
tion. Figure In] depicts the likelihood map of the sky region 
aro und PSR J0218-f 4232 based on the third EGRET cata- 
log ( Hartman ct al. | 1998). The position of the millisecond 
pulsar and the gamma-ray source appear well separated 
from each other. Therefore, the millisecond pulsar PSR 
J0218-I-4232 is very unlikely the counter part of the neigh- 
boring gamma-ray source 3EG J022-t-4253 which is iden- 



tified now with the AGN 3C66A (Hartman et al. 1998) 



7. The isolated millisecond pulsar J2124-3358 

PSR J2 124— 3358 belongs to the small group of iso- 
lated millisecond pulsars. Like PSR J1024— 0719 and 
J1744-1134, it was discovered by |Bailes et al. [ (1997) dur- 
ing the Parkes 436 MHz survey of the southern sky. 
The pulsar has a rotation period of 4.93 ms and a 
proper motion corrected period derivative of P = 1.077 x 
10"^'^ s s"""^, implying an upper limit to the pulsar age of 
P/2P ~ 7.3 X 10^ years and a rotational energy loss of 
E — 3.545 X lO'^'^ erg s~^. Its close distance of about 250 
pc in combination with spin-parameters similar to that of 
PSR J0437-4715 made the pulsar a promising candidate 
for X-ray observations. 

ROSAT HRI observations of PSR J2124-3358 were 
therefore performed on November 3-7, 1995 with an ef- 
fective exposure time of 88 488s. Figure |l2| depicts a 
16 X 16 arcmin^ sub-image of the pulsar field. The pulsar 
is visible as the X-ray point source RX J212443. 5-335841 
at the center of the field. A spatial analysis reveals a po- 
sitional offset between the pulsar's timing position and its 
X-ray counterpart of 4.8", well within the uncertainty of 
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Fig. 11. Confidence regions (50%, 68%, 95%, 99%) of 
the likelihood test statistics for the position of the 
EGRET source neighboring the millisecond pulsar PSR 
J0218-f4232. Also shown is the position of the gamma- 
ray source as listed in the second EGRET catalog (2EG 
J0220-I-4228) and the position of the millisecond pulsar. 
The location of the gamma-ray source is now in agree- 
ment with the position of two AGNs, of which 3C66A is 
foun d to be the like ly counterpart of 3EG J0222-f4253 
(see Hartman et al. 1998 for further details). 



the satellites pointing accuracy. The HRI count rate of 
PSR J2124-3358 is (2.6 ± 0.2) x lO'^ cts/s. 

We applied a timing analysis to 243 photons selected 
from an 8" aperture centered on the pulsar's X-ray posi- 
tion. This cut radius includes approximately 97% of the 



source photons (David et al. 1997). About 10% of the se- 
lected photons are background photons. After correcting 
the data for the satellite's motion and applying a barycen- 
ter correction, a folding of the photon arrival times implies 
strong evidence for the existence of a pulsed signal. The 
resulting light curve is shown in Fig. |l3|. According to the 
H-Test the pulsations are significant at the ~ Aa level (for 
two harmonics). This establishes RX J212443.5-335841 as 
the pulsar's X-ray counterpart and shows the source to be 
the first galactic solitary millisecond X-ray pulsar. 

The pulsar's X-ray light curve shows some indication 
for a double peak as observed in the radio profile at 436 
MHz. Both in the radio and the X-ray profile the phase 
separation between the two peaks appears to be the same. 
Using the bootstrap method the fraction of pulsed X-ray 
photons is estimated to be 33 ± 8%. 

For a comparison between the relative phases of the 
X-ray and radio pulses we have fitted the ROSAT clock 
against UTC for the days 261 - 324 of the year 1995, i.e. 
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Fig. 12. Sub-image of the HRI's field of view based on the 
November 1995 ROSAT observation of PSR J2124-3358. 
The pulsar's X-ray counterpart RX J212443. 5-335841 is 
indicated by an arrow. Several neighboring X-ray sources 
are detected only few arcmin from the pulsar position. 



for a period of 9 weeks. This is long enough to model 
the long term clock drift and to cover the observation 
which was performed during days 307 — 311 of 1995. On 
the other hand we are not limited by clock uncertainties 
which might have happened long before or after the pul- 
sar observation. Calibration measurements are available 
close to the pulsar observation on days 303 and 310 of 
1995. The residuals of our SCC-UTC fit were found to 
be of the order of ±1 ms, corresponding to a la error 
in the relative phase of 20%. In view of this error it is 
not possible to draw a quantitative conclusion concerning 
the phase relation between the X-ray and radio pulse of 
PSR J2124— 3358. However, one may speculate that the 
similarity of both profiles with respect to the phase sepa- 
ration of 0.4 between the pulse components may suggest 
a phase alignment of the dominant radio pulse with the 
weak X-ray pulse. 

Assuming a power-law spectrum with a photon-index 
of a = —2 the detected HRI rate corresponds to an en- 
ergy flux of /j; « 2 — 3 X lO^^'^ erg s^^cm^^ and an X-ray 
luminosity of « 1.5 - 2.2 x 10^°{d/0.25 kpc)^ erg s"^ 
in 0.1—2.4 keV. The latter implies an X-ray efflciency of 
Lx/E w 0.4 — 0.6 X 10^^. The pulsed luminosity is esti- 
mated to be i^p'^ = 6.2 X 10^^ erg s~^. Based on the ra- 
dio dispersion measure of 4.6 pc cm^^ (Bailes et al. 1997) 
and the HI in the Galaxy (Dickey 



Lockman 1990) the 
interstellar absorption is estimated to be in the range 
Nh^2-5x 1020 cm-2. 
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Fig. 13. Light curves of PSR J2124-3358 as observed with 
the ROSAT HRI in 0.1-2.4 keV (top) and the Parkes 
radio telescope at 436 MHz (bottom). Two phase cycles 
are shown for clarity. Both light curves are characterized 
by two peaks which are separated by ^ 0.4 rotations. 



8. Conclusions 

Nine millisecond pulsars are currently detected in the soft 
X-ray domain. Five of them have been identified only by 
their positional coincidence with the radio pulsar and in 
view of their low number of detected counts do not pro- 
vide much more than a rough flux estimate. These objects 
are so faint that the high collecting power of AXAF and 
XMM will be needed to detect enough photons required 
for a detailed spectral and temporal study in the soft and 
the hard band beyond 2 keV. More detailed results are 
found for the other four millisecond pulsars, all of which 
provide interesting empirical information on the pulsar's 
X-ray emission mechanisms. 

PSR 1821-24 in the globular cluster M28 shows X-ray 
pulses up to 10 keV. The sharp peaks in the pulse profile 
and the power-law nature of the spectrum doubtlessly ar- 
gue for a non-thermal origin of the detected emission. The 
alignment between the radio and X-ray pulse component 
adds further support to this interpretation and implies a 
common emission site for the main X-ray and radio pulse 
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components observed at 800 MHz. For PSR J0218+3242 
the sharp peaks found in the X-ray pulse profile also argue 
for a non-thermal origin of the emission. 

The situation for PSR J2124-3358 and J0437-4715 is 
not that clear. For PSR J2124— 3358 there is some indica- 
tion of a double peak structure in the X-ray pulse profile, 
which in terms of pulsed components and pulse phase sep- 
aration implies a similarity between the X-ray and radio 
profile observed at 436 MHz. However, the significance of 
the second peak is not very strong so that further observa- 
tions are required to establish this similarity. Furthermore, 
for broad pulse profiles the shape of a profile itself is not a 
strong indicator for the origin of the detected emission: a 
radiation cone which yields sharp peaks at one aspect an- 
gle may well be seen as hardly modulated away from this 
angle. While sharp peaks indicate non-thermal emission 
processes the reversal - soft modulated emission originates 
from thermal processes - need not be true. The case of 
PSR J0437-4715 is uncertain as well. Here the bandwidth 
limitation of ROSAT and the statistical limitations of the 
ASCA and SAX observations prevent us from conclusively 
discriminating between multi-component thermal spectra 
(thermal polar-cap emission) and a non-thermal origin of 
the radiation. 

Putting the observed emission properties of the de- 
tected millisecond pulsars in a somewhat wider frame 
Becker k. Triimper (1997) found recently that the X-ray 



luminosity of the detected millisecond pulsars show the 
same linear relationship with the same X-ray efficiency as 
the Crab-like pulsars, indicating that the bulk of their 
emission is mainly due to non-thermal processes. The 
occurrence of power-law spectra in PSR B1821— 24 and 
PSR J0437-4715 and the similarity between the radio/X- 
ray pulse profiles seen also for PSR J0218-I-4232 and 
J2124— 3358 may be considered as providing additional 
evidence for a non-thermal origin of the millisecond pul- 
sars' X-ray emission. 

In order to better constrain the origin of the X- 
radiation of the millisecond pulsars reported in this pa- 
per, additional observations with better photon statistics 
are needed. AXAF and XMM, which are designed to have 
higher sensitivity and better spectral, timing and spatial 
performance than ROSAT and ASCA are well suited to 
study millisecond pulsars and to enhance our understand- 
ing of these objects in many aspects. 
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